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It has been suggested that a naked singularity is a good candidate for a strong gravitational 
wave burster. A naked singularity occurs in the generic collapse of an inhomogeneous dust ball. We 
study the odd-parity mode of gravitational waves from a naked singularity of the Lemaitre-Tolman- 
Bondi spacetime. The wave equation for gravitational waves is solved by numerical integration 
using the single null coordinate. It is found that the naked singularity is not a strong source of the 
odd-parity gravitational radiation, although the metric perturbation grows in the central region. 
Therefore, the Cauchy horizon in this spacetime should be marginally stable with respect to odd- 
parity perturbations. 

PACS number(s): 04.20.Dw,04.25.Dm,04.30.Nk 



I. INTRODUCTION 



The study involved in the attempt to detect gravitational waves represents one of the most exciting field in modern 
physics. The detection will provide a decisive test of general relativity and open a new window for astrophysical 
observation. Unfortunately, in spite a great deal of effort, no one has yet succeeded in the detection of the gravita- 
tional waves. However ground-based laser interferometers, such as LIGO, VIRGO, [|| TAMA |^ and GEO600, 
Q are being constructed and will soon enter a stage of data taking. Strong sources of gravitational waves should 
be regions where gravity is general relativistic and where the velocities of bulk motion are near the speed of light. 
Many researchers have theoretically investigated coalescing and colliding black holes and neutron stars as candidates 
for detection. In addition, gravitational waves emitted during stellar collapse to black holes have been studied by 
perturbative calculations of a spherically symmetric background. |^,^ Because of the existence of event horizon, the 
emitted gravitational waves are dominated by the quasi-normal modes of the Schwarzschild black hole in the asymp- 
totic region. Some other realistic candidates exist, e.g., supernovae, rapidly rotating neutron stars, relic gravitational 
waves of the early universe, stochastic gravitational waves from cosmic strings, and so on. 

In this paper we attempt to investigate whether a naked singularity, if such exists, is a strong source of gravitational 
radiation, and, moreover, to understand the dynamics and observational meaning of naked singularity formation. 
Several researchers have shown that the final fate of gravitational collapse is not always a singularity covered by an 
event horizon. For example, in the Lemaitre-Tolman-Bondi (LTB) spacetime, |^,^ a naked shell-focusing singularity 
appears from generic initial data for spherically symmetric configurations of the rest mass density and a specific 
energy of the dust fluid. [^-12| The initial functions in the most general expandable form have been considered, p^ ] 
The matter content in such a spacetime may satisfy even the dominant energy condition. In this case with a small 
disturbance of the spacetime, very short wavelength gravitational waves, which are created in the high density region 
around a singularity, may propagate to the observer outside the dust cloud because of the absence of an event horizon. 
If this is true, extremely high energy phenomena which cannot be realized by any high energy experiment on Earth 
can be observed. In this case information regarding the physics of so-called 'quantum gravity' may be obtained. Also, 
these waves may be so intense that they destroy the Cauchy horizon. In this paper the generation of gravitational 
waves during the collapse of spherical dust ball with small rotational motion is considered. 

Nakamura, Shibata and Nakao |0 have suggested that a naked singularity may emit considerable gravitational 
wave radiation. This was proposed from the estimate of gravitational radiation from a spindle-like naked singularity. 
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They modeled the spindle-hke naked singularity formation in gravitational collapse by the Newtonian prolate dust 
collapse for a sequence of general relativistic, momentarily static initial data. It should be noted that the system they 
considered is different from that considered in this article and that their result is controversial. There are numerical 
analyses that both support and do not support the results of Nakamura, Shibata and Nakao for prolate collapse 
and for cylindrical collapse. JT6| , pT| 

Due to the non-linear nature of the problem, it is difficult to analytically solve the Einstein equation. Therefore, 
numerical methods will provide the final tool. However, its singular behavior makes accurate numerical analysis very 
difficult at some stage. In this article, we investigate odd-parity linear gravitational waves from the collapse of an 
inhomogeneous spherically symmetric dust cloud. Even for the linearized Einstein equation, we must perform the 
numerical integration. However, in contrast to the numerical simulation of the full Einstein equation, high precision is 
guaranteed for the numerical integration of the linearized Einstein equation, even in the region with extremely large 
spacetime curvature. Furthermore, the linear stability of known examples of naked singularity formation is necessary 
as a first step to understand the general dynamics near the "naked singularity formation." 

Recently, Iguchi, Nakao and Harada (INH) studied odd-parity metric perturbations around a naked singularity 
in the LTB spacetime. In INH, it was found that the propagation of odd-parity gravitational waves is not affected 
by collapse of a dust cloud, even if there appears a central naked singularity. When we consider the generation of 
gravitational waves from the dust collapse, we should analyze the perturbations including their matter part. Here 
we investigate the evolution equation of the odd-parity quadrupole mode for metric and matter perturbations. This 
matter perturbation relates to small rotational motion, and it produces gravitational waves during the collapse. We 
follow the derivation of the evolution equations of the perturbations and the numerical method to integrate these 
equations of INH. We investigate the time evolutions of the gauge invariant metric variables at the symmetric center, 
where a naked singularity appears, and at a constant circumferential radius R. We show that the gauge invariant 
variable diverges only at the center and it does not propagate to the outside. 

As we know, the LTB spacetime is one candidate as a counterexample of the cosmic censorship hypothesis (CCH), 
which was introduced by Penrose. ^9 20 1 The CCH is a very helpful assumption, because various theorems on the 
properties of black holes were proved under this assumption. Here the precise formulation and validity of cosmic 
censorship is not our concern. We only consider the situation in which the extremely high-density and large-curvature 
region can be seen by some observer from the gravitational collapse. Such a situation can be regarded as a naked 
singularity in a practical sense since we are not yet able to predict phenomena beyond the Planck scale. However, 
here it should be noted that even from this practical point of view, the stability of the Cauchy horizon with respect 
to gravitational- wave perturbations still has an important physical meaning. Since in the system considered here, the 
mode propagating with the speed of light is the gravitational waves only, the instability of the Cauchy horizon implies 
that the gravitational wave introduces an extremely large spacetime curvature along the null hypersurface near the 
Cauchy horizon from the region near the central singularity. Hence, if the Cauchy horizon in this system is unstable, 
naked singularity formation of this type might be a strong source of gravitational radiation. 

Here we comment on the problem motivated in INH, i.e., "nakedness of the naked singularity." Odd-parity matter 
perturbations are produced by rotational motion of the dust cloud. Therefore the growth of the matter perturbation 
would cause the centrifugal force to dominate the radial motion of the cloud. Hence, an angular momentum bounce 
would occur. If this occurs at the center, the central singularity will disappear. This situation seems to be inevitable 
for the dipole mode of odd-parity matter perturbations. In the case of a spherically symmetric system composed of 
counter-rotating particles, the angular-momentum bounce can actually prevent the naked singularity formation. ]2l| ] 
However, for a system without spherical symmetry, it is still an open question how the non-linear asphericity works 
in the final stage of singularity formation with rotational matter perturbations because of the difficulty of both the 
analytic and numerical treatment. Further, if the initial matter perturbation is sufficiently small, it might still be 
possible that the radius of the spacetime curvature becomes the Planck length. In this case, the behavior of other 
modes of perturbations is crucial to understand the classical dynamics in the region of the Planck scale. As will be 
shown below, since there is no gravitational wave of the dipole mode and the quadrupole mode generated in the dust 
cloud does not propagate to the outside, the Cauchy horizon is marginally stable against odd-parity perturbations 
originating in the aspherical rotational motion of matter. 



The paper is organized as follows: In §0, the basic equations are developed; in §IH, the numerical results are 
presented; these results are discussed in j |lV| ; in we summarize our results. We adopt the geometrized units, in 
which c = G = 1. The signature of the metric tensor and sign convention of the Riemann tensor follow Ref. |22|. 
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II. BASIC EQUATIONS 



We consider the evolution of odd-parity perturbations of the LTB spacetime up to Unear order. The background 
spacetime and perturbation method are described in INH. JlSf In this section we briefly review them. 

Using the synchronous coordinate system, the line element of the background LTB spacetime can be expressed in 
the form 

ds^ = g^^dx^dx'' EE ~dt^ + A^{t, r)dr^ + R^{t, r){d9^ + sin^ 9d(j)^). (2.1) 

The energy-momentum tensor for the dust fluid is 

T^" = p{t,r)u''u'', (2.2) 

where p{t, r) is the rest mass density and is the 4- velocity of the dust fluid. In the synchronous coordinate system, 
the unit vector field normal to the spacelike hypersurfaces is a geodesic, and there is a freedom of which timelike 
geodesic field is adopted as the hypersurface unit normal. Using this freedom, we can always set u'^ — (5q , since the 
4- velocity of the spherically symmetric dust fluid is tangent to an irrotational timelike geodesic field. 

Then the Einstein equations and the equation of motion for the dust fluid reduce to the simple equations 

drR 



(2.3) 

/ N 1 1 dF(r) 



1 1 dj 

8^ R'^drR ~dr 



{mf - ^ - /(r), (2.5) 

where f(r) and F(r) are arbitrary functions of the radial coordinate, r. From Eq. ( ^.4| ), F(r) is related to the 
Misner-Sharp mass function, pSj m{r), of the dust cloud in the manner 



R{t,r) 



F{r) 



m{r) = An / p{t,r)R^dR ^ An / p{t,r)R^drRdr = (2.6) 
Jo Jo 2 

Hence Eq. (|2.5| ) might be regarded as the energy equation per unit mass. This means that the other arbitrary function, 
/(r), is recognized as the specific energy of the dust fluid. The motion of the dust cloud is completely specified by the 
function F{r) (or equivalently, the initial distribution of the rest mass density, p) and the specific energy, f{r). When 
we restrict our calculation to the case that the symmetric center, r = 0, is initially regular, the central shell focusing 



singularity is naked if and only if dfp\r=o < is initially satisfied for the marginally bound collapse, /(r) = 0. 25 
For a collapse that is not marginally bound, there exists a similar condition as an inequality for a value depending on 
the functional forms of F{r) and f {r). ||ll|,||,|5| 

Let us now derive the perturbation equations. The perturbed metric tensor is expressed in the form 

g^Lu = g^Lu + h^^u, (2.7) 

where g^v is a metric tensor of a spherically symmetric background spacetime and is a perturbation. The 
energy-momentum tensor is written in the form 

Tf,. = f^u + 5T^,,. (2.8) 

where T^^ is a background quantity and 5T^i, is a perturbation. By virtue of the spherical symmetry of the background 
spacetime, T^^ can be expressed in the form 

f^.^dx^'dx" = fabdx^dx^ + ^TjfR^{t, r)dn^, (2.9) 

where the subscripts and superscripts a,b, ■ ■ ■ represent t and r, while A, B, - ■ ■ represent 9 and (f>. The odd-parity 
perturbations of /i^^ and (5T^^ are expressed in the form 

/io(i,r)$ri3 

V = I hiit,r)^^s I , (2.10) 



sym h2it,r)xi 



AB 
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io(i,r)$-5 
sym t2{t,r)xT'AB 



(2.11) 



where and xTab odd-parity vector and tensor harmonics associated with the spherical symmetry of the 

background spacetime. |^ We set all the arbitrary constants in the definitions of harmonics to unity. We then 
introduce the gauge-invariant variables defined by Gerlach and Sengupta. [p7[ The metric variables are given by 



ka = ha - T^R^da ( 



The matter variables are given by the combinations 



La — ta T^Lj^ha, 



In the LTB case, the odd-parity gauge-invariant matter variables become 

Lo = p(t, r)U{t, r) and ii = L = 0, 



(2.12) 

(2.13) 
(2.14) 

(2.15) 



where U{t, r) represents the perturbation of the 4- velocity as 5u^ = (0, 0, U{t, 7')$™^). The evolution equation for the 
matter variable (Eq. (3-19) in INH), 



is easily integrated, and we obtain 



1 dJ{r) 



(2.16) 



(2.17) 



^i?2 dr ' 

where J(r) is an arbitrary function depending only on r. From Eqs. (2.4), (2.1?;), and ( 2.17| ), we obtain the relation 



^ C/(r), 

so U{t, r) is independent of the time coordinate t. We introduce a gauge-invariant variable for the metric as 



The metric perturbation variables, kg and fci, are reconstructed from the linearized Einstein equations, 

dr {RS^s) + A{l~l){l + 2)ko^ I&tiAR^Lq, 
dt + ^ (/ - 1) (/ + 2) fci = 0, 



(2.18) 
(2.19) 

(2.20) 

(2.21) 
(2.22) 



by substituting tpg. It was shown in INH that -05 is closely connected to the tetrad components of the magnetic 
part of the Weyl tensor. These components diverge where and only where il^s diverges. From the linearized Einstein 
equations we obtain the linearized evolution equation for the odd-parity perturbation as 



dt 



^dt (i?^.)) - dr (^dr + (^ - 1) (^ + 2) Ats 

( 1 dJ" 



-167r9,. 



V AR^ dr 



(2.23) 



The regularity conditions at the center are also considered in INH. The result for the matter perturbation Lq is 
given by 
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Lo — > icWr'+i + 0(r'+3). (2.24) 

Therefore the matter perturbation variables vanish at the regular center independent of the value of I. The metric 
variable ips behaves near the center as 

i^s — > ^sc{ty~^ + 0{r^) for I > 2, (2.25) 

— > ^sc{t)r + 0{r^) for / = 1. (2.26) 

From the above equations, it is shown that only the quadrupole mode, / = 2, of '0s does not vanish at the regular 
center. 



III. RESULTS 



We numerically solve the wave equation ( 2.23 ) in the case of marginally bound collapse, f{r) — 0, and the quadrupole 
mode, / = 2. We follow numerical techniques employed in INH to integrate the time evolution of the perturbations. 
The numerical code is essentially the same as that in INH, which is tested by comparison with the analytic solution 
for the Minkowski background. 

By virtue of the relation /(r) = 0, we can easily integrate Eq. ( ^.5[ ) and obtain 



R{t,r) 



9F 



1/3 



[^o(r)-^]^/^ 



(3.1) 



where to{r) is an arbitrary function of r. The naked singularity formation time is to = to(0). Using the freedom for 
the scaling of r, we choose i?(0, r) — r. This scaling of r corresponds to the following choice of to{r): 



.3/2 



3VF 



(3.2) 



Here note that, from Eq. (p.3|), the background metric variable, A, is equal to drR- Then, the wave equation (2.23) 
becomes 



1 5Vs 



1 



.drR dlR 



6- 



R 



drR 



f dtR 
dr [ R 



drR 



dips 
dt 



IGtt 



dtdrR \ dtR 1 fdtR 
drR J R 2 \ R 



rd, 



{drR)R^ V {drR)R^ 



(3.3) 



where p(r) = p(0,r) is the density profile at t = 0. We solve this partial differential equation numerically. 

Before getting into the detailed explanation of the numerical techniques, we comment on the behavior of the matter 
perturbation variable Lq around a naked singularity on the slice t — Iq. The regularity conditions of Lq and p 
determine the behavior of U (r) near the center as 



U{r) cx r' 



(3.4) 



This property does not change even if a central singularity appears. However, the r dependence of R and A near the 
center changes at that time. Assuming a rest mass density profile of the form 



p{r) = po + Pnr'' + 



(3.5) 



we obtain the relation 



toir) oc to + tnr"" 



(3.6) 



from Eqs. (2.4) and (3.2), where n is a positive even integer. After substituting this relation into Eq. (3.1), the lowest 
order term is absent from the square brackets of it. Then we obtain the behavior of R and A around the central 
singularity as 



5 



i?(io,r) cx 



(3.7) 



and 



A{to,r) oc ra' 



(3.8) 



on the slice t — to. As a result, we obtain the r dependence of Lq around the center when the naked singularity 
appears as 



Lo{to,r) cx r' 



;-2n+l 



(3.9) 



For example, if I — 2 and n — 2, then Lq is inversely proportional to r and diverges toward the central naked 
singularity. Therefore the source term of the wave equation is expected to have a large magnitude around the naked 
singularity. Thus the metric perturbation variable V's as well as matter variable Lq may diverge toward the naked 
singularity. 

Instead of the {t,r) coordinate system, we introduce a single-null coordinate system, (u,r'), where u is an outgoing 
null coordinate and chosen so that it agrees with t at the symmetric center and we choose r' — r. We perform 
the numerical integration along two characteristic directions. Therefore we use a double null grid in the numerical 
calc ulati on. A detailed explanation of these coordinates is given in INH. By using this new coordinate system, (u, r'), 
Eq. (3.3) is expressed in the form 



du 



a 
R 

a 
2 



2.drR+]^R{dtR)dtdrR-^{dtRfdrR 



d^R 
{drRY 



R 



(1 - dtR) 



?7ra / p{r)U{r) 



dr'^s = 4<^s - 3^^ (1 + dtR) Ips, 



(3.10) 
(3.11) 



wher e the ordinary derivative on the left-hand side of Eq. (3.10) and the partial derivative on the left-hand side of 
Eq. ( 3.11 ) are given by 



— - d + —d - a - " 

du ^ du ^ " 2drR 



dr' 



2drR 



dr, 



dr' = ^^^^dt + dr = {drR)dt 
{dtu)r 



dr, 



respectively. Also, (j)s is defined by Eq. (3.11) and a is given by 

1 



{dtu)r 



(3.12) 
(3.13) 

(3.14) 



We integrate Eq. ( |3.10 ) using the scheme of an explicit first order difference equation. We use the trapezoidal rule. 



sj + l 



■0s 



{{dr'Tps)] + idr'i^s)j+l) , 



(3.15) 



to integrate Eq. (3.11). 

For the boundary condition at the center we demand that t/jg behaves as ipscit) + V's2(i)''^ on a surface oft — const. 
We numerically realize this condition by two-step interpolation. First the values of i/'s are derived at two points on 
the surface oi t — const from the interpolation on the slices of m = const. Next, using these two values, the central 
value of V's is derived from the interpolation on the slice oit = c onst . Another w ay to determine the central value of 
■08 is as follows. We first obtain the central value of 0^ from Eq. (3.10). From Eq. ( ^.11 ) and the boundary conditions, 
the relation of ij^g and <j)s at the center is given by 



— Sdr'R'ips 



(3.16) 



Using this relation the central value of ips is obtained. In our numerical analyses the results of these two methods 
agree well. 
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We assume ■0s vanishes on the initial null hypersurface. Therefore, there exist initial ingoing waves which offset 
the waves produced by the source term on the initial null hypersurface. In INH, it is confirmed that this type of the 
initial ingoing waves propagate through the dust cloud without net amplification even when they pass through the 
cloud just before the appearance of the naked singularity. Therefore those parts of the metric perturbations would 
not diverge at the center. 

We adopt the initial rest mass density profile 



I + exp 



lr_ 

2 r2 



pW=Po „ X , (3.17) 

where po, ri and r2 are positive constants and n is a positive even integer. As a result the dust fluid spreads all over 
the space. However, if r » ri,r2, then p{r) decreases exponentially, so that the dust cloud is divided into the core 
part and the envelope which would be considered as the vacuum region essentially. We define a core radius as 

rcorc=?'l + y. (3.18) 

If we set n = 2, there appears a central naked singularity. This singularity becomes locally or globally naked depending 
on the parameters {po,ri,r2)- However, if the integer n is greater than 2, the final state of the dust cloud is a black 
hole independently of the parameters. Then we consider three different density profiles connected with three types 
of the final state of the dust cloud, globally and locally naked singularities and a black hole. The outgoing null 
coordinate u is chosen so that it agrees with the proper time at the symmetric center. Therefore, even if the black 
hole background is considered, we can analyze the inside of the event horizon. Corresponding parameters are given 
in Table ||. Using this density profile, we numerically calculate the total gravitational mass of the dust cloud M. In 
our calculation we adopt the total mass M as the unit of the variables. 



The source term of Eq. (3.10) 



5(^,0 = (3.19) 



R \ {drR)W- 

is determined by U (r). As m entio ned above, the constraints on the functional form of U{r) are given by the regularity 



condition of Lq- From Eq. (3.4), U{r) should be proportional to r'+^ toward the center. We localize the matter 



perturbation near the center to diminish the effects of the initial ingoing waves. Therefore we define U(r) such that 

r^p(r)U{r) = | (^)' {kf)' < r < (3 2o) 

[ for r > r^, 

where Uq and rb are arbitrary constants. In our numerical calculation we chose rt as rcorc/2. This choice of rb has no 
special meaning, and the results of our numerical calculations are not sensitive to it. 

First we observe the behavior of i/'s at the center. The results are plotted in Fig. |^. The initial oscillations correspond 
to the initial ingoing waves. After these oscillations, tps grows proportional to {to — t)~^ for the naked singularity cases 
near the formation epoch of the naked singularity. For the case of black hole formation, ipg exhibits power-law growth 
in the early part. Later its slope gradually changes but it grows faster than in the case of the naked singularity. For 
the naked singularity cases, the power-law indices S are determined by (Iq — t)'ips/^s locally. The results are shown 
in Fig. H. From this figure we read the final indices as 5/3 for both naked cases. Therefore the metric perturbations 
diverge at the central naked singularity. 



TABLE I. Parameters of initial density profiles, power law indices, and damped oscillation frequencies. 





final state 


po 


ri 


r2 


n 


power index 


damped oscillation frequency 


(a) 


globally naked 


1 X 10"^ 


0.25 


0.5 


2 


5/3 




(b) 


locally naked 


1 X 10"^ 


0.25 


0.5 


2 


5/3 


0.37+0.089i 


(c) 


black hole 


2 X 10"^ 


2 


0.4 


4 




0.37-h0.089i 
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We also observe the wave form of V's along the line of a constant circumferential radius outside the dust cloud. The 
results are shown in Figs. ^-||. Figure ^ displays the wave form of the globally naked case (a), Fig. ^ displays the wave 
form of the locally naked case (b), and Fig. ^ displays the wave form of the black hole case (c). The initial oscillations 
correspond to the initial ingoing waves. In the case of a locally naked singularity and black hole formation, damped 
oscillations dominate the gravitational waves. We read the frequencies and damping rates of these damped oscillations 
from Figs. ^ and|| and give them in terms of complex frequencies as 0.37+0.089« for locally naked and black hole cases. 
These agree well with the fundamental quasi-normal frequency of the quadrupole mode {2Mlo = 0.74734 + 0.17792i) 
of a Schwarzschild black hole given by Chandrasekhar and Detweiler. |28| In the globally naked singularity case (a), 
we did not see this damped oscillation because of the existence of the Cauchy horizon. In all cases the gravitational 
waves generated by matter perturbations are at most quasi-normal modes of a black hole, which is generated outside 
the dust cloud. Therefore intense odd-parity gravitational waves would not be produced by the inhomogeneous dust 
cloud collapse. We should not expect that the central extremely high density region can be observed by this mode of 
gravitational waves. 

We can calculate the radiated power of the gravitational waves and thereby grasp the physical meaning of the 
gauge-invariant quantities, js) To relate the perturbation of the metric to the radiated gravitational power, it is useful 
to specialize to the radiation gauge, in which the tetrad components — ^(0)(0) and ft(6/)((/)) fall off as 0(l/i?), 

and all other tetrad components fall off as 0{1/E?) or faster. Note that, in vacuum at large distance, the spherically 
symmetric background metric is given by the Schwarzschild solution, where hereafter we adopt the Schwarzschild 
coordinates. 



ds" = - 1 - 



2M 



'-If) 



sm 



The relation between the line elements Eq. (2.1) and Eq. ( 3.21| ) is given by the transfer matrix 

1 



dr 



1 - {dtRr 
dR = dtRdt + drRdr. 



[dt + drRdtRdr), 



(3.21) 

(3.22) 
(3.23) 



In this gauge, the metric perturbations in Eq. (2.10) behave as 



/io,/ii = O ( — 



/i2 = w(r - i?,) + 0(1), 



(3.24) 
(3.25) 



where 



i?, =i? + 2A/ln(^^-l 



const. 



Then, the gauge-invariant metric perturbations ( 2.12| ) are calculated as 

fco = -^^«^'^i? + 0(l), 
fci = iu.(i)i? + 0(l). 



(3.26) 

(3.27) 
(3.28) 



where w'^' denotes the first derivative of w with respect to its argument. 

In this radiation gauge, the radiated power P per unit solid angle is given by the formula which was derived by 
Landau and Lifshitz p9] from their stress-energy pseudo-tensor: 



dP _ R^ 
dn ^ 16^ 



I dr 



4 I 9r dT ) 



(3.29) 



For the axisymmetric mode, i.e., m = 0, the above formula is reduced to 



(3.30) 
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2/ + 1 . ^ f(PPi{cost 



where 

^.(^)^^sin^^^^^(^^j . (3.31) 

Then, by using the gauge-invariant quantities and integrating over the all solid angles, the formula for the power of 
gravitational radiation is obtained in the following form: 



where 



Bi = (3.34) 



Using Eq. ( 2.221 ), the radiated power P of the quadrupole mode is given by 



P^^R''[dr{R'i^st. (3.35) 

Figure ^ displays the time evolution of the radiated power P. The radiated power also has a finite value at the Cauchy 
horizon. The total energy radiated by odd-parity quadrupole gravitational waves during the dust collapse should not 
diverge. 

IV. DISCUSSION 

First we con si der t he behavior of the source term S{t, r) around the naked singularity. From the regularity conditions 



and Eqs. (3.4), (3.7), and (B.8), the asymptotic behavior of the source term is obtained as 

S{t,r)(xr^~^ (4.1) 

for t < to and 

5(t,r) cx r'~i"-\ (4.2) 

at t = t(). For example, in the case I — 2 and n — 2, the souce term behaves on t — to as 

S{t,r) (xr~^^/^, (4.3) 

and then it diverges at the center. Thus the divergency of ^ps at the center originates from the source term. To 
confirm this, we numerically integrate the source term along the ingoing null lines with respect to u and estimate the 
central value of i^^. We define this 'estimated' value as 

= / Sit,r)du. (4.4) 



Using Eq. ( 3.16 ) we can define the estimated value of ips as 



We plot it in Fig. ^together with the corresponding ips- The estimated value has the same power-law index of 'ips- 
We conclude that the behavior of ips is determined by the source term in the dust cloud. 

We next consider the stability of the Cauchy horizon. We found that the metric perturbation produced by the 
source term does not propagate outside the dust cloud, except for quasi-normal ringing. The source term, which 
controls ips, does not diverge at the Cauchy horizon. Therefore ips should not diverge at the Cauchy horizon and 
should not destroy it. Then, even if odd-parity perturbations are considered, it will not be the case that the LTB 



9 



spacetime loses its character as a counterexample to CCH due to Cauchy horizon instabihty. Also, it does not seem 
that such collapse is a strong source of gravitational waves. 

In this paper, we have dealt with the marginally bound case. For the case of non-marginally bound collapse, the 
condition of the appearance of the central naked singularity is slightly different from that in the above case [^4|j25[ 
and hence there is the possibility that the behavior of ips m this case is different from that in the marginally bound 
case. However, it is well known that the limiting behavior of the metric as i — > to{^) is common for all the cases: |29| ] 



We conjecture that the results of the perturbed analysis for the non-marginal collapse would be similar to the results 
for the marginal bound. 



We have studied the behavior of the odd-parity perturbation in the LTB spacetime including the matter perturba- 
tion. For the quadrupole mode, where gravitational waves exist, we have numerically investigated the wave equation. 
For the case of naked singularity formation, the gauge- invariant metric variable, diverges according to a power law 
with power index 5/3 at the center. This power index is closely related to the behavior of the matter perturbation 
around the center. We have also observed ips at a constant circumferential radius. For the globally naked case, we 
cannot see intense gravitational waves propagated from the center just before the crossing of the Cauchy horizon. For 
the locally naked case, we have confirmed that there exist quasi- normal oscillations. As a result, we conclude that 
the type of singularity changes due to the odd-parity perturbation because Tps diverges at the center. However, the 
Cauchy horizon is marginally stable against odd-parity perturbations. 

At the final stage of the collapse, the effects of the rotational motion are important and the centrifugal force 
might dominate the radial motion. If this is true, the central singularity would disappear when an odd-parity matter 
perturbation is introduced. For the dipole mode, such a situation seems to be inevitable. However, we should note 
that it is a non-trivial and open question how non-linear asphericity affects the final fate of the singularity-formation 
process. Further, in the case of initially sufficiently small aspherical perturbations, the radius of spacetime curvature 
at the center might reach the Planck length, and hence there is still the possibility that the naked singularity is formed 
there in a practical sense. However, as our present analysis has revealed, since the Cauchy horizon is stable with 
respect to odd-parity linear perturbations, there is little possibility that this collapse is a strong source of odd-parity 
gravitational waves. 

There remains important related works to be completed. The first problem is to consider the even-parity mode in 
which the metric and matter perturbations are essentially coupled with each other. We are now investigating this 
problem. Finally, we should consider the non-linear effects to complete this analysis. This problem will be analyzed 
in the future. 
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FIG. 1. Plots of V's at the center as a function of the time coordinate t. The solid line represents the globally naked case 
(a), the dashed line represents the locally naked case (b), and the dotted line represents the black hole case (c). 
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FIG. 2. Plots of the local power indices {to — t)tps/ip3- The solid line corresponds to the globally naked case (a), and the 
dashed line corresponds to the locally naked case (b). Both of them approach a value near 5/3. 



01 




-IxlO-' 



220 240 260 280 300 320 
(-) t 



-lxl0-» 



-2x10-8 



-3x10-8 



322 



324 



326 



328 



(b) 



FIG. 3. Plots of 'ips for the globally naked case (a) at R = 100. In (a) the left hand side oscillation originates from the initial 
ingoing wave. In (b) we magnify the the right-hand edge, which is just before the Cauchy horizon. The dotted lines represent 
the time at which the observer at i? = 100 intersects the Cauchy horizon, which is determined by numerical integration of the 
null geodesic equation from the naked singularity. 
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FIG. 4. Plots of for the locally naked case (b) at R = 100. In (a) the left hand side oscillation originates from the initial 
ingoing wave. After this oscillation, the damped oscillation dominates, and this part of the wave form is magnified in (b). 
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FIG. 5. Plots of for the black hole case (c) at 7i = 100. In (a) the left-hand side oscillation originates from the initial 
ingoing wave. After this oscillation, the damped oscillation dominate, as depicted in (b). 
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FIG. 6. Plots of the radiated power P for globally naked case (a) at J? = 100. The horizontal axis is the out-going null 
coordinate u. At the Cauchy horizon, this coordinate has the value uq. 
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